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ABSTRACT 

The contribution from the residuals of the foreground can have a significant impact 
on the temperature maps of the Cosmic Microwave Background (CMB). Mostly, the 
focus has been on the galactic plane, when foreground cleaning has taken place. How- 
ever, in this paper, we will investigate the possible foreground contamination, from 
sources outside the galactic plane in the CMB maps. 

We will analyze the correlation between the Faraday rotation map and the CMB tem- 
perature map. The Faraday rotation map is dependent on the galactic magnetic field, 
as well as the thermal electron density, and both may contribute to the CMB tem- 
perature. We find that the standard deviation for the mean cross correlation deviate 
from that of simulations at the 99.9% level. Additionally, a comparison between the 
CMB temperature extrema and the extremum points of the Faraday rotation is also 
performed, showing a general overlap between the two. Also we find that the CMB 
Cold Spot is located at an area of strong negative cross correlation, meaning that it 
may be explained by a galactic origin. 

Further, we investigate nearby supernova remnants in the galaxy, traced by the galactic 
radio loops. These super nova remnants are located at high and low galactic latitude, 
and thus well outside the galactic plane. We find some correlation between the Fara- 
day Rotation and the CMB temperature, at select radio loops. This indicate, that the 
galactic foregrounds may affect the CMB, at high galactic latitudes 
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1 INTRODUCTION 

The study of the cross-correlations between the Faraday ro- 
tation and the Cosmic Mic rowave Background (C MB) tem- 
perature, was started in jPineen fc Colej |2004| ). and has 
showed itself to be of remarkable importance in the field 
of CMB science. 

The first, third, fifth and seventh year data re- 
leases from the Wilkinso n Mic row ave Anisotropy Prob e 
( ^WMAP) dB ennett et aU '2003"), (Hinshaw etH] 120071 ). 
IjHinshaw et al . 2009), (Jarqsik, ct al. 2011) provide the 
most precise measurements of the CMB temperature dis- 
tribution. The precision of these measurements will be sig- 
nificantly improved with the coming data from the ongoing 
PLANCK mission. In order to be used for precision cosmol- 
ogy, the CMB data must be cleaned from the residuals of 
the galactic foreground. The focus for this cleaning has been 
the galactic plane, as this includes the most powerful con- 
tributions from the galaxy. However, in this article, we will 
investigate the possible foreground contamination at high 



and low galactic latitude. In particular, we will look at the 
effect that the interstellar medium (ISM), as well as nearby 
supernova remnants, might have on the CMB. 
The turbulent ISM consist of hot {Tg > W'^K) low den- 
sity (rie ~ 10~^cm~^) ionized gas, filling about 70 — 80% 
of the Galactic volume and interspersed by cold neutral 
and relatively dense clouds with a high density of HI atoms 
( < 10~^cm~^). These HI regions are the places for intensive 
star formation, which can turn the corresponding HI clouds 
into the short lived HII zones. The corresponding transi- 
tion was observed in the blue compact galaxy Henize 2-10, 
by detection of a compact ( < 8pc) 1 mjy radio sources 
dKobulnickv fc Johnsonl|l999l ). If similar HI-HII transitions 
occurs in the Milky Way, these clouds could manifest them- 
selves in the CMB range of frequencies as a feature of the 
optical depth for the CMB photons, creating negative spots 
in the directions of the clouds, due to the SunyaevZeldovich 
effect. 

While the influence of the ISM clouds on the CMB signal 
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could have profound implication for the CMB-cosmology 
and the ISM-science, this aspect of the problem attract very 
little attention in the literature. There is global mass bal- 
ance by evaporation, ionization and condensation between 
supernova energy i njection and radiative coo li ng in the ISM 
(see fo r instance , l|McKee fc Ostrikeij [TqTtI I , l|Melioli et all 
I2OO9I ). l|Stilll2009l ') and corresponding references therein). 
Additionally, the galactic sky contains several radio spurs, 
which can be joined together into loops. These are called 
the galactic radio loops, and since the 70s, six loops 
have been identified, commonly referred to, as loop I- 
VI (for a review see (|Borkall2007l )). Loop I-IV are desig- 
nated as the mai n radio loops, and were d is covered first 
llLaree et al.lll962ll . (jQuiglev fc Haslamlll965l ). l|Large et. all 
1 19661). whereas loops V and VI, were only confirmed later 



i Milogradov-Turin fc Smithlll973[ ). The origin for the radio 
loops axe the quasi- spherical remnants of very old, nearby 
supernov ae explosions, a nd the ensuing shock fronts (see for 
instance (|Urosevidl2002l )V The shock front contains an en- 
velope of HI atoms, and violent mixing of warm and hot 
gasses takes place inside the shock area. The subsequent syn- 
chrotron emission from these supernova remnants (SNR), 
at the galactic radio loops may contribute to the CMB- 
radiation in form of the foreground residuals. We will inves- 
tigate this, by analyzing correlations between the magnetic 
field strength and the CMB temperature, in the loop areas. 
The goal of our paper is to show that the HI-HII clouds can 
interfere with the CMB signal, and to investigate a possible 
connection between the galactic radio loops and the CMB 
temperature. To test these ideas we looked at the cross- 
correlations between a WMAP temperature patch of the 
sky r(r) and the Fa raday rotation depth $ (r), using re- 
cently available data (|Oppermann et aLll201ll ). 
The outline of the paper is the following. First we introduce 
the concept of Faraday depth, as well as our cross correla- 
tion coeffecient, for circular areas on the sky. Then we intro- 
duce the two datasets that we use, the data for the Faraday 
depth, and for the CMB, as well as our masking procedure. 
In the next section we construct the map of correlations, and 
present areas of extreme correlations, as well as investigate 
connections between Faraday depth and CMB temperature 
extremes. We hereafter compare the correlation maps to sim- 
ulations, analyzing the significance of our findings. We then 
turn to the specific area around the Cold Spot, and inves- 
tigate the significance of the correlations in the area, using 
various tests. Next, we start on the analysis of the loop areas 
in the sky, and look at the correlations there between CMB 
temperature and Faraday depth, as well as between polar- 
ized data and Faraday depth. Finally, in the last section, we 
summmarize the paper, and draw our conclusions. 



2 CMB TEMPERATURE-FARADAY 

ROTATION DEPTH CROSS-CORRELATION. 

The first test which we would like to implement, is the com- 
parison of the CMB signal and the depth of the Faraday 
rotation "1> for various areas of the sky. Historically, the cross- 
correlations of the CMB and the F araday rotation map was 
detected in l|Dineen fc Coles! |2004| ) for selected patches of 
the sky, mainly due to incompleteness of the last one. In 
a recent paper (jOppermann et ahl I2OIII ) , the most exten- 



sive catalog of Faraday rotation data of compact extragalac- 
tic polarized radio so urces from the VLA sky survey NVSS 
(|Condon et al. )ll998l). has been assem bled. 
The authors ( Qppermann et aLll201lf ) introduce the concept 
of Faraday depth, which depends on position and is indepen- 
dent of any astrophysical source. The Faraday depth, corre- 
sponding to a position at a distance ro from an observer, is 
given by a line of sight integral, 



$(r-o) = 



27rm?c4 



drnE{r)Be{r) — 



(1) 



over the thermal electron density and the line of sight 
component of the magnetic field Br. Here, e and are the 
electron charge and mass, c is the speed of light, and 8 and 
<j) are the polar and azimuthal angles of the polar system of 
coordinates. From Eq[T]one can see a remarkable feature of 
the function I {6, 0). If the magnetic field is nearly constant 
within some coherence lenght Xb, the Faraday depth $(S, (f)) 
and I {6, (p) are proportional to the Thompson optical depth 

In order to look for correlations between the CMB tempera- 
ture measurements and Faraday rotation measure, we define 
the cross-correlation map rs{l,b,R) as follows 



rs{l,b,R)^ 



(2) 



where e = 10" , is used to avoid the possible divergence in 
the calculation. C is the circular region around {l,b) with 
radius being R. x and y are the mean values of Xi and yi 
in the kq75pl masked maps. In our analysis, Xi denotes the 
Faraday map and yi denotes the CMB temperature map. It 
is clear that —1 < < 1, and > indicates the positive 
correlation while Vs < shows a negative correlation. If R — 
360°, rs{l, b, R) becomes independent of {l,b), and reduces to 
the well-defined full-sky correlation coefficient. When R <^ 
360°, rs{l,b,R) describes the local correlation of Faraday 
rotation map and CMB map around the location (Z, b) with 
scale R. 



3 DATA USED IN THE ANALYSIS 

3.1 Galactic Faraday rotation map 

Faraday rotation measurements of the extragalactic radio 
sources provide the tracers of the Galactic magnetic field. 
When the polarized radiation propagates through a mag- 
netized plasma, the Faraday depth is given by Eq. [1] In 
the recent work (|Oppermann et al]|201ll ) , the authors recon- 
structed the full-sky Galactic Faraday depth maps by using 
the point source Faraday depth measurements of the NRAO 
VLA SKY Survey (NVSS) catalog as well as other surveys. 
In the reconstruction, in order to divide out the most obvi- 
ous large scale anisotropy introduced by the presence of the 
Galactic disk, the authors defined the dimensionless signal 
s{l, b) as 



sil,b) 



b) 
p{b) ' 



(3) 



where <;/>(/, b) is the potential Faraday depth map. The vari- 
ance profile p{b) is a function of Galactic latitude only. 
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We show the reconstructed signal map s{l,b) and the de- 
rived Faraday depth map (j>(l,b) in Fig. [1] and Fig. [2] which 
have the HEALPix l|Gorski et al.l 120051 ) resolution parame- 
ter JVaide = 128. In both figures, we see a large region with 
quite large uncertainties in the southern sky below the dec- 
lination of —40°, which is caused by the lack of data in this 
regio n due to the position of the observing telescope VLA 

i see JOppermann et al.ll201ll 'l. (|Oppermann et al.ll201ll l and 
Condon et al. 1 1998l 'l for the details). In our following anal- 
ysis, we shall mask this region to reduce the uncertainties. 



3.2 Data of CMB temperature fluctuations 

The WMAP instrument is composed of 10 differencing as- 
semblies (DAs) span ning five frequencies from 23 to 94GHz 
l|Bennett et al.|[2003l ). The internal linear combination (ILC) 
method has been u sed by the WMAP team to generate th e 
WMAP ILC maps | Hinshaw et "ani2007l ), (|Gold et aLlbOllT ). 
The 7- year ILC (ILC7) map is a weighted combination from 
all five original frequency bands, which are smoothed to a 
common resolution of one degree. In this paper, we shall 
consider the ILC7 map for the analysis. 
Also, we will analyze the 7-year WMAP foreground reduced 
Q-band (Q7), V-band (V7) and W-band (W7) maps, as weU 
as the co-added maps of these frequency bands. The linearly 
co-added maps are constructed using an inverse weight of 
pixel-noise variance. We adopt two maps with different com- 
binations of frequency bands: V-band and W-band (written 
as 'VW7') and Q-band, V-band and W-band (written as 
'QVW7'). Note that all these WMAP data have the same 
resolution parameter A^sidc = 512. In order to cross-correlate 
with Faraday depth maps, we degrade them to a lower res- 




Figure 3. The kq75pl mask. 



olution A'sido ~ 128 for the analysis. 

To exclude the effect of various contaminations, in the anal- 
ysis we apply the KQ75pl mask, which is a combination of 
KQ75y7 mask given by WMAP team and a mask which re- 
jects the regions with large uncertainty in the Faraday signal 
map, i.e. as{l,b) > 0.7 (see the right panel in Fig. [T]). The 
mask is shown in Fig. [3l 



4 REGIONS OF STRONGEST CORRELATION 

We now construct the rs{l, b, R) map by the following steps: 

1) We degrade the CMB temperature anisotropy maps 
(ILC7, Q7, V7, W7, QVW7, VW7) to the low resolution 
A^sidc = 128. 

2) For the given Faraday map (0-map or s-map) and the 
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degraded CMB map, we mask them by kq75pl mask. 

3) Based on these two masked maps, we calculate the values 
of X and y. 

4) Choosing an R value, and using Eq.((2|, we construct 
rs{l, b, R) map. 

Note that for every C, we only take into account the 
unmasked pixels. In the output rs{l,b,R) map, we set 
ra{li,bi, R) = if more than 10% pixels are masked in the 
C around {k, bi). 

In Figs. |4] and [5] (upper panels), we plot the rs{l, b, R) maps 
with R = 5° and R = 7° where the Faraday (/!>-map and the 
ILC7 CMB map are used. Also, we calculated the same cor- 
relation for the full sky map (masked by the KQ75pl mask), 
resulting in an overall correlation of 0.0504. 
In order to investigate the Faraday-CMB cross correlation 
further, we define two statistics: the mean value fs and the 
standard deviation , which are calculated from the un- 
masked rs{l,b,R) maps. Specifically, we calculate the stan- 
dard deviation as 

= ^^E(-»--)^ (4) 

where Xi is the i'th pixel value out of A'^ in the region, and 
X is the average pixel value for the unmasked pixels. 
Now, for each rs(l,b,R) map, we calculate the mean value 
fs and standard deviation , and only show the areas with 
strongest correlation, i.e. \rs — 7^i| > 2(7,,^. The results are 
presented in Figs. [4] and [5] (bottom panels). We have also 
checked that we obtain similar results, if we replace the 
Fararay (^-map with the s-map, or replace the ILC7 map 
with the other CMB maps, or consider other cases with 
R = 3° and 9°. 

According to Fig.s|4]and[5]we have several zones of pe- 
culiar cross-correlations between the Faraday rotation map 
and the CMB temperature map. It is interesting to see, 
whether or not the extrema of the Faraday depth coincide 
with extreme temperatures in the CMB map. To investigate 
this, we present in Figs. [B] [3 |S] and [5] the CMB temperature 
isolines overlaid upon the Faraday depth map for 4 such ar- 
eas. Clearly, we see several areas, where the concentration 
of isolines and the Faraday depth extrema overlap. This in- 
dicate a connection between the galactic magnetic field and 
the extreme temperatures of the CMB. 



4.1 Comparing with simulations 

In order to study the global properties of the cross- 
correlation between Faraday rotation maps and the WMAP 
data, we shall compare the rs maps based on WMAP data 
with those of the random CMB simulations. 
The cosmology in the random Gaussian simulation is based 
on the ACDM with th e cosmological param eters from 
WMAP 7-year best-fit l|Komatsu et all boill ): Qbh^ = 
0.02255, n^h'^ = 0.1126, = 0.725, h = 0.702, r = 0.088, 
A|j = 2.430 X lO"'', Us = 0.968. First, we simulate the 
CMB maps for each frequency channel by considering the 
WMAP beam resolution and instrument noise for each chan- 
nel, and then co-add them with inverse weight of the full-sky 
averaged pixel-noise variance in each frequency to get the 
simulated QVW7 and VW7 maps. To simulate the ILC7 
maps, we consider the 1° smoothing resolution and ignore 




Figure 4. Top: rs{l, 6, R) with R = 5° based on Faraday t/i-map 
and ILC7 map; Bottom: same as top one, but here we have set 
= if |rs - I < 2ar^ . 




-0.715 ^^^^^^^^^^^^^^^^^^^^^^^^^ t^^^^^^^^^^^^^B 



Figure 5. Top: rs{l, b, R) with R = 7° based on Faraday (/i-map 
and ILC7 map; Bottom: same as top one, but here we have set 
rs = if |rs — r"s I < 2a-r, . 
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Figure 6. Top. Isolines of the ILC 7 CMB temperature fluctua- 
tions on the Faraday rotation depth map, for the zone with Galac- 
tic coordinates I = 37°, b = 54°. The isolines start at —0.025 mK, 
and go down with a step of —0.025 mK. Bottom. The same, as 
the top, but for all levels from O.OmK and up, with a step 0.0125 
mK. Note the non-symmetric color scale, in which magenta indi- 
cate negative signal, while, starting from blue, the Faraday depth 
is positive. 



Figure 7. Similar to Fig. [6] but at Galactic coordinates I = 
125°, fc = 47° 



the noises. In the simulation, we assumed that the temper- 
ature fluctuations and instrument noise follows the Gaus- 
sian distribution, and do not consider any efi'ect due to the 
residual foreground contaminations. We have analyzed one 
thousand WMAP simulated data sets in the same way as 
the real data set. 

Since the physical mechanism for the cross-correlation be- 
tween Faraday rotation maps and the CMB temperature 
anisotropy maps is quite complicated, which depends on 
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Figure 8. Similar to Fig. |6] but at Galactic coordinates I = 
44°, b = 14° 



both the magnetic distribution and the thermal electron dis- 
tribution, we cannot expect to know how the location and 
amount of highly correlated areas in the WMAP deviate 
from those of random Gaussian simulations. 
We compare the ?^^ and dr^ from the WMAP data with 
those of 1000 random simulations, and calculate the follow- 
ing probability functions 

Q = Prob(lr,"l > |r,|), /3 = Prob(cr^^ > a,J, (5) 
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Figure 9. Similar to Fig. [6] but at Galactic coordinates I = 
200°, 6 = 54° 



where r^" and a^^ are the values for the WMAP data, while 
7^^ and (Jr^ are those for 1000 CMB realizations. If there is no 
correlation between Faraday rotation maps and the WMAP 
data, we expect that the values of a and /3 should be close 
to 0.5. A larger deviation from 0.5 indicates a stronger cor- 
relation. 

We calculate a and P values for all the combinations of 
the Faraday maps (i/i-map, s-map) and CMB maps (ILC7, 
QVW7, VW7, W7, V7, W7), by considering the cases with 
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Figure 10. The values of r~s for ILC7 and Faraday t/i-map (red 
line) compared with those for 1000 simulations. Upper left is for 
R = 3° , upper right is for R = 5° , lower left is for R = 7° , and 
lower right is for R = 9° . 




Figure 11. Same as Fig. llOl but the ILC7 map is replaced by the 
W7 map. 



J? = 3°, 5°, 7°, 9°. We find that for the statistic r^, the 
deviations from the random simulations are quite small, i.e. 
a € [0.66, 0.74] is satisfied for all cases. These are clearly 
shown in Figs. [TU] and 1111 

However, we find that the values a,.^ of WMAP data are 
quite smaller than those of simulations, i.e. the /3 parame- 
ters are dramatically small for all the cases (see Tables fOI 
and 14. II and Figs. I12land ll3|l . For instance, for all the cases 
related to ILC7 CMB map and Q7 map we always have 
P < 0.001. So we conclude that the cross-correlations for 
these WMAP data deviate from those of random realiza- 
tions at 99.9% confidence level. From these Tables, we also 
find that the deviation from random realizations becomes 
smaller for the higher frequency bands. For the cases related 
to CMB W7 data, /3 becomes ~ 0.05 for R = 3° and ~ 0.005 
for R = 9° . Even so, we also find that the cross-correlations 
for the WMAP data deviate from those of random realiza- 
tions at more than 94.5% confidence level. The dependence 
on the CMB frequency band can be explained as follows: in 
the higher frequency bands, the contaminations of Galactic 
magnetic fields and/or the thermal electrons (such as the 
residual synchrotron emissions) are smaller than those in 
the lower frequency bands. 



Table 1. Results of <Tr, and P, at R = 3°, 5°, 7° and 9°, for 
Faraday t/i-map and various CMB data. 





ILC7 


QVW7 


VW7 


Q7 


V7 


W7 


< (3°) 
/9 (3°) 


0.438 
0.000 


0.372 
0.008 


0.362 
0.020 


0.379 
0.001 


0.360 
0.004 


0.357 
0.050 


< (5°) 
/3 (5°) 


0.349 
0.000 


0.294 
0.002 


0.287 
0.003 


0.299 
0.000 


0.285 
0.001 


0.283 
0.016 


< (7°) 
13 (7°) 


0.294 
0.000 


0.246 
0.001 


0.240 
0.002 


0.251 
0.000 


0.239 
0.001 


0.237 
0.005 


< (9°) 
13 (9°) 


0.255 
0.000 


0.213 
0.001 


0.208 
0.002 


0.217 
0.001 


0.206 
0.001 


0.205 
0.005 


Table 2. 

Faraday s- 


Results of (7r^ and 13, at R 
map and various CMB data. 


= 3°, 5° 


,7° and 


9°, for 




ILC7 


QVW7 


VW7 


Q7 


V7 


W7 


< (3°) 
13 (3°) 


0.435 
0.000 


0.369 
0.009 


0.359 
0.020 


0.375 
0.000 


0.358 
0.006 


0.354 
0.054 


< (5°) 
13 (5°) 


0.344 
0.000 


0.290 
0.001 


0.282 
0.003 


0.295 
0.000 


0.281 
0.001 


0.278 
0.010 


< (7°) 
13 (7°) 


0.289 
0.000 


0.243 
0.001 


0.236 
0.001 


0.247 
0.000 


0.235 
0.001 


0.233 
0.004 


< (9°) 
/9 (9°) 


0.250 
0.000 


0.209 
0.001 


0.203 
0.001 


0.213 
0.001 


0.202 
0.001 


0.201 
0.006 



5 CORRELATIONS AT THE COLD SPOT 

It is most interesting, that from Figs. |4] and [5] we find a 
strong negative correl ation around the C MB Cold Spot at 
{I = 209°, b = -57°) (jBennett et al.ll201ll ). 
Originally, the Cold Spot ( CS) w as detected in a 
wavelet analysis JVielva et al.l |2004| ). (|Cruz et all |2005| ). 
JCavon et al.l bOOsF of the first-year data release from 
WMAP. It is apparently inconsistent with the assump- 
tion of statistically homogeneous Gaussia n fiuctuations, an d 
its existence has later been co nfirmed (ICruz et al.l |2006| ) , 
(|Cruz et all |2007| ). Previously, (|Rudnickl ^o'dtF have ana- 
lyzed the NVSS maps and have discovered a decrease in the 
space density of radio sources in the CS region. An indepen- 




Figure 12. The values of ar^ for ILC7 (red line) comparing with 
those for 1000 simulations. Upper left is for R = 3°, upper right is 
for R = 5° , lower left is for R = 7° , and lower right is for ij = 9° . 
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Figure 13. Same with Fig. 1121 but ILC7 map is replaced by W7 
map. 




Figure 15. Same as Fig. [H but for Ts{l,b,R) based on (/)-map 
and CMB W7 map. 



Figure 14. Histogram for tiic (/>-map and ILC7 based rs{l, b, R) 
maps witli R = 3° (upper left), R = 5° (upper right), R = 7° 
(lower left), R = 9° (lower right). In all panels, the red lines 
indicate the values around the CMB Cold Spot. 



resolution, we define the mosaic cross-correlation K{1, b) as: 



dent study showed that the CS, m ay be a simple statistica l 
deviation due to systematic effects (|Smith fc Hutererll2010l ) , 
or peculiaritie s of the low multipol e tail of the CMB map 
at 2 < / < 10 ijNaselskv et al.l 120101 '). There are a few mod- 
els of the pri mordial (c osmolosical) origin o f the CS, based 
on te xtures (|Cruz et al.. ,2008. '), large voids (|Das fc Spergell 
I2009D etc. 

In order to quantify our findings regarding the CS, we inves- 
tigate what the rs-values is for the CS area, compared to the 
rest of the sky. In Figs. 1 141 and 1 151 we present histograms of 
the masked ra{l, b, R) maps for four different radii (i? = 3°, 
7? = 5°, i? = 7° and i? = 9°) for the ILC and "W7 cases. 
The red line in the histograms corresponds to the value of 
rs around the CS. All histograms show that at the CS 
is smaller than mean value for the map at more than 2(t 
level. Clearly, the region around the CS show an anoma- 
lously strong anticorrelation, compared with the rest of the 
map. 

We now turn to the mosai c correlation method 
llVerkhodanov fc Khabibullin"all2010l ). to study the correla- 
tion of the Faraday depth signal in the CS zone with the 
WMAP ILC data. 

For any two maps with the same pixelization and angular 



K(l,b) = 




- < X 




" i6n(i,b) 




2 




»^ 




ien{i,6) 




2 


5Z < 2/» 


>r 



(6) 



ien{;,6) 



where l,b denote the coordinates of the zone ^l{l,b) = 
S° X S° , with pixel numbers i G SI and < .. > denote the 
average over all pixels i in the zone b). 
Using K{1, b) from Eq|6l we now produce a mosaic correla- 
tion of the ILC map and the Faraday depth measure map. 
In Fig. [16] we present a cut of 30° x 30° , around the cold 
spot area, using a correlation window of Sl(/, &) = 2° x 2°. 
Here, the blue spots corresponds to anticorrelational zones 
of the ILC and Faraday depth signal, and the red spots in- 
dicate zones with positive correlations. We calculated statis- 
tics of these correlational coefficients and compared it with 
1000 simulations of the Faraday rotation and CMB model 
correlational data. The histogram on Fig. [17] show the dis- 
tribution of the coefficients in the CS region (left), and for 
the total sky (right). Notice the two-peaks distribution, dif- 
ferent from the expected for an uncorrelated signal. Such a 
distribution is connected with the close positions of maxima 
and minima in the ILC and Faraday depth maps. 

Finally, similar to the procedure above, we test the area 
around the CS for overlap between the extrema of the Fara- 
day depth and of the CMB temperature. The results are 
presented in Fig. 1181 It is clear to see, in the top panel, that 
the concentration of isolines and the areas of high and low 
Faraday depth overlap well. The high temperature isolines in 
the right panel does not fit as good, although there is some 
overlap in the lower right area of the sky patch. Nonethe- 
less, these results are even stronger, than those from Figs. [6] 
through [9] and underline the peculiarity of the high corre- 
lations in the CS region. 



6 CORRELATION AT THE GALACTIC 
RADIO LOOPS 

We now turn to areas on the map, at high and low galac- 
tic latitude, covered by the galactic radio loops. The shock 




Figure 16. Mosaic correlation map of the ILC signal with the 
Faraday rotation depth in the CS region. The window is 2° X 
2°. (left panel). Right panel. Smoothed map of the correlation 
between the ILC signal and the Faraday rotation depth. 




Figure 17. Histogram of the correlation coefficients in the CS 
region (the left panel), and for the entire CMB sky (the right 
panel). The correlation window is 2° X 2°. The black solid line is 
the normalized number of pixels for the ILC map. The gray area 
demonstrate the ±la-scatter produced for 1000 simulations. 



front of the super nova remnants (SNR), outlined by the 
galactic loops, affects the observed galactic magnetic field. 
The magnetic fields at the shock front may be the origin of 
synchrotron radiation, possibly affecting the CMB. To in- 
vestigate, whether or not the area of the SNRs affect the 
CMB, we cross correlate the Faraday <?!>-map and the ILC7 
map, at the loop locations. The areas in the sky occupied 
by the loops are shown in fig. [19] with the kq75pl mask 
as background, and in table [SI with the limiting longitudes 
and latitudes taken from (jBorkaj|2007|) . The area covered by 
loop I, is overlapping the well known northern galactic spur, 
which is very clearly visible at 45 MHz in Fig. I2UI The data 
for the 45 MHz map is taken from (|Guzmanll201ll ). 

In Fig. [21] we overlay the areas of the radio loops with 
the Ts maps from Figs. [4] and [5] It is interesting, that several 
of the areas of high and low correlation are located inside the 
radio loops. Most interesting are loop II A (in the middle 
lower part) and loop VI A & B (right upper part), which 
coinside very well with areas of negative correlation. Of areas 
with positive correlation, loop V B and somewhat C (the 
second and first small loop, left of II A respectively) stand 
out. It is clear, that there is no preference for either negative 
or positive correlation in the loops, as will also be shown 
below. However the loops do, for the most part, contain 
areas of strong correlation. 

We now investigate the correlation between the Faraday 
(;/>-map, and the CMB ILC7 map, and the V-band map in 
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Figure 18. Similar to Fig. \E\ but for the Cold Spot region {I = 
209°, b = -57°). 




Figure 19. The 6 galactic radio loops, with the kq75pl mask as 
background. Dark red: Loop 1 A & B. Yellow; Loop II A & B. 
Light green: Loop III A, B & C. Cyan: Loop IV A. Blue: Loop V 
A, B & C. Dark blue: Loop VI A & B. 



particular, in the areas of the galactic radio loops. Again we 
use the method of local cross-correlation, introduced earlier, 
only for square areas on the sky this time. In order to mini- 
mize the contribution to the correlation, caused by noise in 
the two maps, we smooth them, before calculating the cross 
correlation. We select smoothing values of 1°, 5° and 10°, for 
both the (j> map, as well as the ILC7 and V-band map. We 
then apply the kq75pl mask to the </>-map and the ILC7 and 
V-band map, which we have scaled to the same resolution as 
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Table 3. The galactic longitude (Z) and latitude (b) in degrees, 
for the loops I- VI. 



Loop 




I n 




Loop I 


A 

B 


[360,327] 


ri Q 7Qi 
[10, (fS\ 

[67,78] 


Loop II 


B 


[57,30] 
[195,130] 


r '^n ml 

[ OU, J-U 1 

[-70,-2] 


Loop III 


A 
B 
C 


[180,135] 
[135,110] 
[110,70] 


[2,50] 

[40,55] 

[6,50] 


Loop IV 


A 


[325,285] 


[55,72] 


Loop V 


A 
B 
C 


[189,178] 
[147,133] 
[90,80] 


[-25,-13] 
[-50,-39] 
[-39,-24] 


Loop VI 


A 
B 


[215,205] 
[207,196] 


[29,40] 
[6,32] 




Figure 20. The sky at 45 MHz, with the north galactic spur 
clearly visible. Color scale is chosen to emphasize the spur struc- 
ture. The blue circular area to the left of the map is not covered 
in the dataset. 



the mask. We only consider pixels, not covered by the mask, 
for the calculation of the cross correlations. Note, that one 
of the loop V areas (A) are almost completely covered by 
the mask, resulting in extreme values of the correlations. 
Finally, we compare the cross correlation value for each loop 
area, with the cross correlation value of 100 simulated CMB 
maps, and count the amount of simulations having a greater 
value of rs. The entire procedure de scribed above, was 
cond ucted using the GLESP program (|Doroshkevich et al.l 
[iooi), with a resolution of fix = 201, meaning that the map 
had 201 pixels in the ^-direction, and 402 in the 6'-direction. 
The results from the cross correlations can be seen in tables 
[H] and [51 for smoothings of 1°, 5° and 10°. Also included in 
the tables, is the number, S, of simulations with a higher 
value of the cross-correlations. 

From the table of results, it is immediately obvious, 
that the smoothing have a big effect on certain loops. 
Because of the contribution from areas outside the loops 
during the smoothing process, the small area loops are 
naturally most affected (see for instance loop III B). 
Further, there is not much difference between the ILC7, 




Figure 21. The contours of the 6 galactic radio loops, overlaid 
with the Ts map for i? = 5° (top) and R = 7° (bottom). 



Table 4. The cross correlation (r^) for the galactic radio loops, 
and the amount of simulations with a higher cross correlation 
(S), for three values of smoothing. This is compared to 100 sim- 
ulations. 



Loop 




Smooth 1° 


Smooth 


5° 


Smooth 10° 






rs 


S 




S 




S 




A 


0.534 


15 


0.737 


8 


0.807 


5 


Loop I 


B 


0.760 


5 


0.888 


10 


0.961 


17 




A 


-0.323 


93 


-0.409 


89 


-0.456 


85 


Loop II 


B 


0.335 


9 


0.377 


10 


0.401 


14 




A 


0.283 


7 


0.445 


7 


0.573 


6 


Loop III 


B 


0.108 


37 


0.553 


19 


0.903 


5 




C 


0.127 


35 


0.328 


29 


0.615 


21 


Loop IV 


A 


0.330 


33 


0.392 


43 


0.447 


50 




A 


0.999 


46 


-0.675 


54 


-0.996 


54 


Loop V 


B 


0.068 


48 


0.600 


28 


0.930 


1 




C 


0.208 


51 


0.259 


51 


0.165 


54 




A 


-0.0311 


49 


0.119 


49 


0.649 


40 


Loop VI 


B 


0.049 


50 


0.233 


46 


0.692 


27 



and the V-band, meaning, that the ILC method does not 
cause any of the correlation effects. We also see, that for 
almost all loops, we have an increase in the r^, when we 
increase the smoothing degree, as the contribution from 
noise diminishes. However, only some of the areas show a 
systematic increase in significance compared to the random 
simulations. Of particular interest are the areas covered by 
loop I A, and loop III A, as these areas cover a relatively 
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Table 5. Similar to tableO but for the V-band only. 



Loop Smooth 1° Smooth 5° Smooth 10° 







rs 


S 




S 




S 


Loop I 


A 
B 


0.532 
0.731 


15 

5 


0.754 
0.887 


8 

10 


0.820 
0.963 


3 

17 


Loop 11 


A 
B 


-0.308 
0.327 


93 
9 


-0.402 
0.367 


88 
12 


-0.459 
0.384 


85 
15 


Loop 111 


A 
B 
C 


0.276 
0.110 
0.117 


8 

37 
35 


0.447 
0.580 
0.347 


7 

19 
29 


0.574 
0.920 
0.601 


6 
1 

22 


Loop IV 


A 


0.345 


31 


0.482 


36 


0.578 


46 


Loop V 


A 
B 
C 


0.999 
0.162 
0.219 


44 
43 
51 


-0.027 

0.669 

0.239 


53 
20 
52 


-0.990 

0.930 

0.0798 


53 
1 

54 


Loop VI 


A 
B 


-0.053 
0.026 


49 
51 


0.092 
0.231 


49 
46 


0.655 
0.746 


40 
23 



large area of the sky, and consistently exhibit a strong 
correlation compared to simulations. Further, as Fig. I2UI 
shows, the boundary of Loop I could be extended left, 
to more accurately overlap the north galactic spur. From 
FigED we see, that this extension would incorporate two 
bright red spots into the correlation coefficient, and thus 
increase the level of cross-correlations above the threshold 
mentioned above. 

Although the results are not very strong, this nonetheless 
shows, that the magnetic field in the powerful radio loops 
may affect the CMB photons. 

Further, we have also investigated the correlations in the 
loop for the CMB polarization. We have used the E-mode 
polarization from the K-band for the test, and comp ared to 
simulations created via Healpix (jCorski et al.|[2005l ). Apart 
from the different maps, the procedure for the correlations 
in the loops is exactly the same, as described above. 
In Fig. [22] we show the E-modes from the K-band. Notice 
the spur-like structure, extending from the left of the 
galactic centre and extending northwards. This have a 
remarkable resemblance to a similar structure, visible in 
roughly the same area of the map in Fig. [1] though it is 
not clearly visible in the (ji-map (Fig And addition- 
ally, the area coincide reasonably well with the area of 
Loop I, as seen in Fig. 1191 making it particularly interesting. 

The results of the correlation test for the polarization 
are presented in table (6) In overview, many results seem 
significant at 1° smoothing, but becomes non-significant at 
greater levels of smoothings. Also interesting is the seem- 
ingly high levels of correlation, where several areas are dis- 
tinctly anti-correlated. 

The area mentioned above, at Loop I (A), like several other 
loop areas, shows a significant anti-correlation at 1° smooth- 
ing, but at greater smoothings the correlation value quickly 
drops to insignificant levels. This could be explained by the 
fact, that the spur-like structure in Fig. [22] seems to be rel- 
atively thin. At greater levels of smoothing, the feature be- 
comes indistinct, and thus the significance drops. 
A more detailed analysis of these correlations will be the aim 




-0.33 ^^^^^^^^^^^^^^^^^^^K +1.11 



Figure 22. The map of E-modes in the K-band. Notice the spur- 
like structure, extending from the left of the galactic centre and 
extending northwards. 



Table 6. Similar to table [6] but for the K-band, E-mode polar- 
ization, compared with simulated polarization data. 



Loop 




Smooth 1° 

Ta S 


Smooth 5° 

rs S 


Smooth 10° 

Ta S 


Loop 1 


A 
B 


-0.712 
-0.622 


99 
90 


-0.732 
-0.763 


83 
77 


-0.726 
-0.831 


78 
73 


Loop 11 


A 
B 


0.099 
0.143 


36 
28 


0.120 
0.150 


43 
40 


0.154 
0.151 


44 
43 


Loop III 


A 
B 
C 


-0.036 
-0.584 
-0.603 


56 
98 
100 


-0.031 
-0.649 
-0.687 


49 
78 
85 


-0.052 
-0.708 
-0.770 


49 
76 
84 


Loop IV 


A 


-0.629 


98 


-0.759 


84 


-0.857 


78 


Loop V 


A 
B 
C 


-0.076 

0.526 

0.529 


58 

3 

14 


-0.068 

0.615 

0.545 


58 
22 
37 


-0.066 
0.677 
0.544 


56 
27 
38 


Loop VI 


A 
B 


-0.642 
-0.590 


89 
91 


-0.648 
-0.622 


64 
73 


-0.643 
-0.643 


63 
68 



for future work. Presumably, when higher quality polariza- 
tion maps, as well as better data for the Faraday depth will 
be available, the results may improve. 



7 CONCLUSION 

In this paper, we have investigated possible foreground con- 
taminations of the CMB, at high and low galactic latitudes. 
We have cross correlated CMB maps, with the Faraday Ro- 
tation map, to identify areas of high or low correlation. 
Specifically, in our first test, we selected circular areas, of 
various radii for our comparison test. The idea is, that the 
Faraday Rotation is affected by thermal electrons, as well as 
the strength of the galactic magnetic field. A large value for 
the Faraday Rotation, at a given area, can imply that we 
have a high column density of thermal electrons, and thus 
create a cold area in the CMB. Our cross correlation test 
for the ILC map showed that there was a considerable anti- 
correlation between the Faraday rotation map and the ILC 
map, at the Cold Spot (CS), as well as in a number of other 
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areas, all well outside the galactic plane. The correlations in 
the CS area were compared to correlations over the rest of 
the sky, where the results showed the CS area as highly an- 
ticorrelated. We further tested the area around the CS with 
the mosaic correlation method, with a 2° x 2° window. 
In a histogram of the the number of pixels contra correlation 
coefficients for the mosaic correlation, we found two signif- 
icant peaks compared to random simulations, at negative 
correlation values, in good agreement with the results for 
the cross correlations. 

Additionally, we compared the Faraday depth and the CMB 
temperature isolines, in an area around the CS and at 4 other 
locations of high correlation value. We found in general, that 
many of the CMB temperature extrema, and Faraday depth 
extrema overlapped, substantiating the theory that the mag- 
netic field affects the CMB temperature. Note that this were 
all for areas well outside the galactic plane. 
Lastly, we compared our maps of correlation with simulated 
maps, and found that the standard deviation for the mean 
cross correlation in the ILC case, deviated from random sim- 
ulations at the 99.9% level. 

For our second test, we repeated the cross correlations, but 
this time at very specific areas of the map. We focused on the 
galactic radio loops, which is the remnants of supernova ex- 
plosions. The shock fronts of the supernovas, outlined by the 
galactic radio loops, contains strong magnetic fields, possi- 
bly giving rise to synchrotron radiation. Thus we cross corre- 
lated the Faraday Rotation maps, in the areas of the galactic 
radio loops, and compared these with random simulations. 
During the correlations, we smoothed both maps to various 
degrees, to avoid contamination from the noise. The results 
for these tests were not as strong, as for the unspecified area 
test above, but nonetheless interesting. For the ILC case, we 
found the most powerful results for the heavy smoothing of 
10 °, where loop 1(A) and loop III(B) both had a significance 
at the 95% level. The area occupied by loop III(A) was also 
interesting, having a relatively high significance, irrespective 
of the smoothing. We repeated the test for the V-band, but 
found similar results, showing that the process of creating 
the ILC map, is not the cause of the correlations. Finally, we 
also investigated the correlation between the Faraday depth 
map and the K-band E-mode polarization, in the loop areas. 
The results showed, that we had significant ant icorrelat ions 
for the case of 1° smoothing, but insignificant correlations 
for greater levels of smoothing. 

Comon to all the tests are, that at current, the Faraday maps 
used are not of the highest quality, and thus the results for 
the comparisons may improve in the future. 
In conclusion, we find that foreground contamination at high 
galactic lattitudes, should be taken very seriously, and if not 
properly filtered, may cause anomalies in the CMB. 
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